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We investigate proposals of how the form factor approach to compute cor- 
relation functions at zero temperature can be extended to finite temperature. 
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interacting theories require a more severe change of the form factor program. 
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1 Introduction 



The computation of correlation functions is one of the central objectives in quantum 
field theories. In general, this can only be achieved by means of perturbation theory 
in the coupling constant. Nonetheless, in 1+1 space-time dimensions many exact 
results are known, in particular at zero temperature. At present, one of the most 
successful approaches in this direction is the form factor program. Originally this 
method was developed to compute correlation functions for massive models at zero 
temperature |1|, . More recently, it has also been demonstrated, that the approach 
can be employed successfully for the computation of massless correlation functions 
H for vanishing temperature. However, in a realistic set-up of a physical experiment 
one needs to know such functions in the finite temperature regime. Computing for 
instance physical quantities from linear response theory one requires the response 
function in form of the canonical two-point correlation function at finite temperature 

(AO(x,t)AO'(x',t , )) T , (1.1) 

where AO(x, t) = 0(x, t) — (0(x, t)) T . For the static case, e.g., for the computation 
of electric and magnetic susceptibilities this goes back to [|[]. Taking for instance 
the local operators 0(x,t) and 0'(x,t) to be the current J(x,t), the dynamical 
response is needed to compute the conductivity by means of the celebrated Kubo 
formula ||. 

In U a proposal was made to adapt the form factor approach to finite temper- 
ature computations. In there it was demonstrated for some operators of the Ising 
model that the modified approach indeed reproduces the expected results for the 
temperature dependent correlation functions, even when boundaries are included. 
In order to make the method meaningful several technical assumptions were re- 
quired to eliminate various infinities. Since not all of them can be justified in an 
entirely rigorous fashion, further evidence is desirable to support the working of 
the proposed prescription, even for the Ising model. The proposal is very appeal- 
ing, since apart from the introduction of a dressing function, the main structures 
of the (T=0)-form factor approach are perpetuated. It was argued in M for the 
Ising model, that the dressing functions admit an interpretation in terms of density 
distribution functions. This observation was taken up in [0] and it was conjectured 
that the interpretation should also hold for interacting theories. Some checks which 
support the validity of this conjecture for the one-point function were presented in 
0. Shortly afterwards, doubts were raised in || on the working of these formulae 
for the two-point function, albeit only a counter example which required a chem- 
ical potential was provided. An additional controversy arose thereafter about the 
nature of the "dressing function" (see equation ( [2.15 )) which has to be employed in 



the context of the one-point function. In 0, it was proposed to employ the on-shell 
free energies, rather than the pseudo-energies obtainable from the thermodynamic 
Bethe ansatz (TBA) as suggested in f7|. For the one-point function evidence was 
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provided in [[H]] that the proposal in |9j appears to be incorrect. No explicit claims 
concerning the two-point function were made in ||. 

The main purpose of this manuscript is to contribute to this debate and provide 
further evidence for the (non)-validity of the various proposals. In regard to the 
importance of ( |1 . 1| ) , we want to focus especially on the study of the expressions for 
the two-point functions. So far the only few explicit computations using the dressed 
form factor approach may be found in || || . 

Our manuscript is organized as follows: In section 2 we outline the various 
proposals made so far to evaluate temperature dependent correlation functions in 
the massive as well as in the massless regime. We investigate these proposals for 
two models: the complex free Fermion/Federbush model (section 3) and the scal- 
ing Yang-Lee model (section 4). In section 5 we state our conclusions. In the 
appendix we assemble various properties of functions which occur throughout our 
computations. 



2 Temperature dependent correlation functions 

We start by providing a concise review of the main features of the prescription to 
compute temperature dependent n-point correlation functions by means of form 
factors. In general a temperature state is described by a density matrix p and 
the expectation value for observables is the trace over the product of this matrix 
with the observables. Taking to be eigenstates of the Hamiltonian H with 
eigenvalues E^, and as usual (3 = 1/kT with k being Boltzmann's constant and 
T the absolute temperature, the temperature dependent n-point function for the 
observables 0\ . . . O n is defined as 

(0 1 (x 1 ,t 1 )...O n (x n ,t n )) T := I^V^ (i/;\0 1 (x 1 ,t 1 )...O n (x n ,t n )\i/;) . (2.1) 

As usual, this expression is normalized by dividing with the partition function 

Z = Xy^(V#> , (2-2) 

i> 

which ensures that the trace over the density matrix p = e~@ H /Z is one. The central 
idea of the "dressed" form factor program is now, as in the finite temperature case 
n, |2| , to reduce the computation of the expansion ( |2.1| ) to the computation of form 
factors 

*y = (0(O)|V>) • (2.3) 
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This is achieved simply by the insertion of (n — 1) complete states = 1. 

Suppressing for compactness the space-time dependence of the operators, this reads 

(0 1 ...0„> r =i Yl e-^o^ \0 1 \^ 1 }^ 1 \0 2 \^ 2 }...(^ n _ 1 \O n \%}. 

(2.4) 

Thereafter one needs a meaningful prescription to relate matrix elements of the 
form (ip'\ O \ip) to those where the vacuum is on the left, as in ( |2.3| ), and a shift 
operator u tjX 0(x,t) = f(x,t)O(0) which moves the operator to the origin. In the 
following we will take x M = (— ir, 0), which implies the restriction r < (3 in order to 
ensure that exp(— (/3 + r)H) is a trace class operator. 

Let us now specify the states to be multi-particle states of the form 

m = \zl{e l )zl{e 2 )...zl n {e n )), (2.5) 

where the operators Z^{9) are creation operators for a particle of type H as a 
function of the rapidity 9. These operators are assumed to satisfy the Faddeev- 
Zamolodchikov algebra |TJ Z\{9 X )Z]{9 2 ) = Sf j {9 l2 )zl{9 2 )Z\{9 1 ) with S being the 
scattering matrix depending on the rapidity difference 8i 2 = 9i — 9 2 . The prescrip- 
tion O \if>) -> (0\tpifj f ) then reads 

(ZM . . . Z» n (d'n)\ 0(x,t)\Zl(9 n ) . . . Zl{0i)) = 

(0(x,t)\Zl n (6 n ) . . . Z\ x {Q x )Z-^) . . . Z ftn (6'-)) (2.6) 

all contractions 

with 9~ = 9 — m + ie and e is an infinitesimal quantity and ft being the anti- 
particle of fi. After shifting the operator to the origin, the r.h.s. of ( |2.6[ ) involves 
the n-particle form factors, which we denote as 

FS^»*~{p x , . . . A) = (O(0)| ZliO,) . . . Zl(9 n )) . (2.7) 

These form factors have to satisfy various properties |TL |2|, such as Watson's equa- 
tions 

Fn^ """ M,Ml+1 """(. . . ; Q i: . . .) = F,f' "' M * +l/v "(- • • , 9 i+1 , Oi, . . .)S llilJi . +1 {6i ji+ x) , (2.8) 
F^-^(9 1 + 27ri,...,6 n ) = 1 ° i F^-^(9 2 ,...,9 n ,9 1 ) , (2.9) 
Lorentz invariance 

F^-^(9 U ...,9 n ) = e ^ F o\^..^ dl + X ,...,9 n + \), (2.10) 
and the so-called kinematic residue equation 



ResF^f *"*»{0 + iir, 9 , 9 X . ..9 n )= i[l - ^\{ S^(9 0l )]F°^-^(9 v ...,9 n ). 

(2.11) 



i=i 
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In ( [2.1 0|) s is the Lorentz spin of O and A G C an arbitrary shift. In equations ( |2.9| ) 



and ( [2. 11|) 7^ is the factor of local commutativity defined through the equal time 



exchange relation of the local operator O(x) and the field O^y) associated to the 
particle creation operators ZUO), i.e., O l _ l (x)0(y) = 7^ 0(y) 0^(x) for x 1 > y l . 
It is the singularity lim e ^o ZJ l (9)Zp / (9~) — > 00, implicit in ( [2. lip , which is the rea- 
son for the presence of the e in (|2.6| ). The renormalization prescription which 
eliminates these divergencies is outlined in |T^] (see also P, [71). Other types of sin- 
gularities arise from the contractions of terms like Z^(9)Zl(9)Z u (9')Zl(9') -> 5 2 {0). 
As demonstrated explicitly in || such terms are absorbed in the (re)-normalization 
factor Z. In this way the one-point function (TJ, ||, [7], ||, |||, [Kj 



n=0 Pi-Hn 

x (Z^Oi) . . . Z, n (9 n )\ 0{r)\Zl n {9 n ) . . . Z^)) r < I, (2.12) 

becomes a meaningful expression. In many cases this function is an important 
normalization factor, however, for the reasons mentioned in the introduction we 
will focus our attention on the two-point function. It results to || 

n 



n=l flx-Hn 1=1 

xF^-^(9 u ...,9 n ) [F°>^(9 h . . .,9 n )]* r < ^.(2.13) 



This formula requires several explanations and comments: We dropped here as 
usual another infinity coming from n = in the infinite sum. The sum over the /i 
extends over particles and holes. This is understood in the way that to each particle 
type present at zero temperature one may associate a hole, such that each term at 
zero temperature which is summed over n-particles is mapped into 2 n -terms at finite 
tempearture. According to ||, the form factors involving holes may be constructed 
from the ones of particles by an Z7r-shift 

F 0| Ml ...particle... / , n ^ i . . . ? 0. ? . . . n ) = F 0|M 1 ...hole... Mn ^ i ^ q, _ ^ _ ^ ^ 

The origin of the hole interpretation is thus the crossing of particles from bra to ket 
by means of (|2.6| ) together with the renormalization prescription. As explained in 
detail in , in comparison with the one-point function there are additional singular 
terms emerging in the two-point functions. In the expansion they occur in terms in 
which the product of the two form factors associated to O and O' involve a different 
amount of particles. These terms are just dropped, which could be a possible source 
for the difficulties we encounter below in interacting theories. 
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The functions fi(9,T) are the so-called filling fractions 



fi(0,T) 



1 -£«(()) exp[- ei (0,r)] 



(2.15) 



involving the functions Ei(9,T). For the non-interacting case, treated in ||, this 
function is taken to be the on-shell energy divided by the temperature Ei(9,T) = 
rrii/T cosh 9. When extending the validity of these formulae to the interacting case, 
one may speculate on the nature of Ei(9,T). In it was proposed to interpret 
Ei(9,T) as the pseudo-energies, which may be determined by means of the thermo- 
dynamic Bethe ansatz equation |L3] 



E,{ 



rm. 



coshfl + ^.b^m/^) 



(2.16) 



As common we denote here f = m/T, mi —>■ mi/m, with m being the mass of 
the lightest particle in the model. By (/ * g) (9):= 1/(2tt) / d9'f(9 - 9')g(9') we 
denote the convolution of two functions and <Pij{9) = — id In Sij(9)/d9. In contrast, 
when specializing the statement in [|7| to fermionic statistics, i.e., Su(0) = — 1 
, it was suggested therein to take instead the free on-shell energies divided by 
the temperature also in the interacting case. More generally, this means that the 



eigenvalue E$ of the Hamiltonian in 



is either taken to be the free on-shell or 



the pseudo-energy. In each case, the pseudo-energies of the holes are simply the 
negative of the ones of the particles 



£ho\e(9, T) = —E 



particle I 



?,T). 



(2.17) 



A simple but key property satisfied by the two-point function is the Kubo-Martin- 
Schwinger (KMS)-condition |14j]. Assuming to have a time evolution operator uo t at 
disposal, it is easily obtained from the trace properties of the temperature dependent 
correlation function 



(u t OO') T = (0'Lu t+if3 0) T (0{x, t)0'(x', t')) T = (0'{x', t')0{x, t + i(3)) 



T 



(2.18) 

For more detailed discussion on this formula see e.g., JT^]. For the choice = 
(— ir, 0), as in ( |2.13| ), this condition reads 

(u t OO') T = (0'cu t _pO) T & (O(r)O'(0)) T = (O'(O)O(r - /3)) T . (2.19) 



Noting that (O{r)O'{0)) T = (£>(0)£>'(-r)) T , it is clear that fl27Lj) indeed satisfies 
the condition (|2.19|) , provided that the form factors obey 

r ~i * 

riO In X holes m X particles / /i n "\ rpO' \nxholes mxparticles//! n \ 

r n (pi, ■ ■ ■ , V n +m) r n v (Pi . . . , V n+m ) = 

T "I * 

r-iO|nxparticles mxholes//) n \ zT'O'lnxparticles mxholes//i n \ /o on\ 

r n {Ui,...,U n+m ) t n ' {V 1} . . . , 9 n+m ) . (2.2UJ 
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However, from one of Watson's equations ( |2.9| ) and Lorentz invariance ( |2.10| ) we 
easily derive that in general one picks up a factor exp iix(s' — s) on the r.h.s. of ( [2.20|) , 
where s, s' are the Lorentz spins of O and O' , respectively. Consequently, there is 
no problem with KMS when (s' — s) G 2Z like for instance when the two operators 
coincide. The latter case was treated in |], |7|, §, |9], Despite the fact that the 
KMS condition puts some structural restrictions on the two-point functions, it is 
not constraining enough in what more precise functional details concern, e.g., it 
can not shed any light on the controversy [0, |^, 10] about the precise nature of the 
dressing function ( |2.15|) . However, it dictates the two functions £ A1 (^,T) appearing 
explicitly in (|2.13|) and ( [2.15| ) to be identical. 

To clarify further the working of (|2.13|) it would be highly desirable to compute 
this functions for more explicit models. Unfortunately there are not many tem- 
perature dependent two-point functions known from alternative approaches which 
one could compare with in order to settle the issue. Nonetheless, various limits are 
known which one can take as benchmarks. 



2.1 The massless limit, conformal correlation functions 

In conformal field theory various methods have been developed to compute cor- 
relation functions. At vanishing temperature the most celebrated approach is the 
one which exploits the structure of the Virasoro algebra, such that the correlation 
functions obey certain differential equations |l(J. To include the temperature is also 
fairly simple in this case. It is achieved just by mapping the observables from the 
plane to the cylinder, z -> exp(2yrT^), O(z) -> (2irT)- A °e- 27T ™ Ao O($) 



= (27rT) A 0+ A , + A 0+ A , (O( Zl )O'(z 2 )) T=0 (0{z x )0'{z 2 )) T=Q . (2.21) 



T 



By construction Q2.21]) satisfies the KMS condition, provided (O(zi)O'(z 2 )) T=0 = 
(O' (z 2 )0(zi)) T=Q . Alternatively one can get some further information on this func- 
tion by exploiting the KMS condition on one of the holomorphic sectors by adopting 



a proposal made in [18]. Ignoring normal ordering one obtains 



OC 



(0(& 1 )0'(<& 2 )) T = E exp[-27rzT(iM + iM0][O„,C>fc] • (2.22) 



-oo k= 



For this to hold we only need to assume that for and 0'(ii} 2 ) exist Fourier- 

Laurent mode expansions of the form 

Eoo 
esxp(-2imiT'd)O n . (2.23) 
n=— oo 

In addition one makes use of the fact that the time evolution is governed by u t O = 

e itL Q e -itL ^ with 

Lq being the zero mode generator of the Virasoro algebra. If 
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then furthermore the commutator [Lq, O n ] = —nO n holds (this is true for instance 
for O n = L n the modes of the energy-momentum tensor, O n = J% the modes of a 
Kac-Moody current, O n = <f) n the modes of a primary field, O n = G n the modes 
of an N = 1 supersymmetric fielcf]), the relation (|2.22|) is derived immediately. 
To make contact with ( |2.21[ ) one needs of course to incorporate a proper normal 
ordering prescription. 




<00' > 

dressed form factors 



^"particle ^' ^"holc ^ 



T = 



filling fractions, TBA 



m = 



< OO' > 

m.T-0 



standard form factors 



F F 

L V L R 



F F 

RR RL 

F F 

L L> L R 



m = 



OO' 



massless dressed form factors 



T = 



OO' 



map to cylinder 



massless form factors 




Figure 1: Two-point correlation functions in various limits. 

Alternatively, we may compute the correlation functions by using the (dressed) 
form factors related to the massless theory. The prescription of taking the massless 
limit was originally introduced in Jl9| within the context of a scattering theory. It 



consists of replacing in every rapidity dependent expression 9 — > 9 ± a, where an 
additional auxiliary parameter a has been introduced. Thereafter one should take 
the limit o — > oo, m — > while keeping the quantity m = m/2exp(<r) finite. For 
instance, carrying out this prescription for the momentum yields p± = ±mexp(±#), 
such that one may view the model as splitted into its two chiral sectors and one can 
speak naturally of left (L) and right (R) movers. Hence, having a function depending 
on the rapidities of n particles, it will be mapped into 2 n related functions 



lim f(6 1 + x i( j, ...,6 n + x n a) = /„!,...,„„ (0 ls • • • , 0„), R , 

a—>oo V i — It IOr Hi — -p 



V; = L for 



x, ; = - 



(2.24) 



*In fll7| it was shown, however, that supersymmetry and temperature seem to be incompatible 
concepts. Only the vacuum states admit the implementation of supersymmetry. 
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For the scattering matrix this means every massive amplitude is duplicated S(9) = 
Sll{9) = S RR {9) and in addition one obtains the amplitudes connecting the two 
chiral sectors S R l/lr(9) = lim^oo S(9 ± 2a). 

In Q this prescription was also carried out for expressions of form factors. In that 
case each n-particle form factor is turned into 2™ n-particle form factors. Note that 
when considering (|2.24 ) for form factors this does in general not lead to the same 



expressions as when taking the scattering matrices Sll(0), S rr (9), Srl(8), Slr(0) 
and determining the form factors thereafter in the usual fashion. This two ways of 
carrying out the limit only commute for form factors associated to operators whose 
Lorentz spin is vanishing. 

In order to be able to formulate the analogue of the expression fl2.13| ) for the 



massless case one also requires the massless version of the dressing function, i.e., the 
massless analogues of the pseudo-energies. For this one can use once more the above 
recipe, such that the TBA-equations ( |2.16| ) are replaced by the same equations with 
S(9) -> S LL (9),S RR (9),S RL (9),S LR {9) and f m; cosh 9 -> f m;exp6>. The working 
of this was confirmed in the analysis of JTDJ. 

Having now outlined the prescription to compute the massless (temperature de- 
pendent) correlation functions, a non-trivial check is constituted by the comparison 
with ( |2.21| ). A more constricting check is to start with ( |2.13|) and carry out the 



massless and zero temperature limit. That is checking the commutativity of the 
diagram in figure 1. In particular the massless limit can possibly shed more light 
on the issue of different dressing functions. 

3 The complex free Fermion/Federbush model 

We will now demonstrate the working of the previous approach with various ex- 
amples. Let us start with the complex free Fermion case (the complex version of 
the Ising model) with Virasoro central charge c = 2 x 1/2. One reason for not 
considering directly the simpler case of self-conjugate Fermions is that all formulae 



presented in this section also hold for the more general Federbush model |22], |23[| . 
This means we can also regard the results of this section to hold for more exotic 
statistical interaction. In general, the free Fermion is particularly attractive to start 
with, since for many operators the higher n-particle form factors are vanishing such 
that the infinite series in ( |2.13| ) terminates. For the self-conjugate Fermion several 
two-point correlation functions at zero temperature have been computed, e.g., in 



|20| , pTf . Despite the simplicity of the model, only for very few operators massive 
temperature dependent correlation functions have been evaluated by means of form 
factors 0. In order to put formula ( j2.13|) on firmer ground it is therefore desirable 



to check its working first of all for a wider range of operators. 

Taking the operator O and O' to be the current, we will present all four cases 
illustrated in figure 1 in some detail. The current-current correlation functions are 
particularly interesting, since they occur explicitly in the application within the 



S 



Kubo formula. Adopting the notation of our recent exposition |23] , we consider the 
correlation function involving one of the chiral currents J = J° ± J 1 . The only 
non- vanishing massive form factors for these operators are 



—mme 



-j-i+t 

+ 2 



(3.1) 



In the following we shall focus on the mutual correlator of J = J~ . Out of the four 
cases depicted in figure 1, the conformal case at zero temperature is the easiest to 
treat and hence a good starting point. According to the massless limit prescription 
( ggg ), we obtain 



r 2,RR\ 

j?J\ii ( 
r 2,LL\ 

rpJ\ii 
1.LL 



(0,9) 



r 2,RR\ 



-2nimexp(6 + 6)/2 



7 J\ti , 

2,RL\ 

7 J\n , 
2,RL\ 



0. 



(3.2) 

(3.3) 
(3.4) 



such that (|2.13j ) yields 



(J(r)J(O)) 



m=0,T=0 



dOdO 



(2tt; 



2 exp[-rm(e" + e y )}(e b 



(3.5) 



Recalling that the current has scaling dimension Aj = 1, Aj = 0, this agrees 
of course with the leading order term of the well-known conformal ?7(l)-current- 
current two-point correlation function k/r 2 for level k = 1, see e.g., [Q. Note, that 
as it should be, the auxiliary parameter rh has vanished in the final expressions. 
Raising now the temperature we can use the same expressions for the form factors, 



but the proposal (|2.13|) dictates that we have to dress them with the massless version 
of the filling fractions 



f±(0,T) 



1 



1 + exp(=prh/Te e 



(3.6) 



for particles (+) and holes (— ). The values are in agreement with ( |2.17|) . With this 
we compute 



(J(r)J(O)) 



m=0,T 



Am 2 TX 2 y j 

7T 2 T 2 

sin 2 (titT) 



dQdQ 

(2^' 



r )U9,T)f v (6,T)e 



-r rh(fj,e s +ue s ) 



(3.7) 



The result ( ^.7[ ) can of course also be obtained directly from the mapping ( |2.21| ) and 
the correlation function at zero temperature (|3.5|). Making now the model massive, 
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we employ instead of (|3.2|) and ( |3.3| ) the form factors ( |3.1| ) and evaluate 



(J(r)J(O)) 



m,T=0 



2 2 
m 7T 



(27T 



2 exp[— rm (cosh + cosh#)](e + ) 



m 2 [lf 1 (rm)]'' 



(3. 



where i^i(x) is a modified Bessel function (see appendix). Using the limiting be- 
haviour ( |A.2|) , we recover as expected the conformal correlation function ( |3.5|) in 
the limit m — > 0. Considering now the massive finite temperature regime, we have 
to include in the previous computation the massive dressing function 



f±(0,T) 



1 + exp(=Fm/T cosh 6) 
Then we compute according to (|2.13|) 



(3.9) 



(J(r)J(O)) 



2 2 

m it 



(27T) 



hS )/„(0, T)f v (9, r) e - rm ^ cosh0+i/cosh ^ 



m 2 LK"/" (m, r, T) 



n 2 



(3.10) 



The functions K^(m, r, T) are defined in the appendix. The restriction on the argu- 
ments of the Bessel functions (|A.1|) reflects the conditions r < 1/T in ( 2.13|) . With 
the help of ( |A.10| ) we recover the expression ( |3.8|) for T — > 0. Taking instead first 
the limit m — > in (|3.10| ), we reproduce with ( |A.5 ) the previously computed con- 
formal correlator ( |3.7| ). In conclusion this means the different methods to compute 
(J(r)J(O)) for several mass and temperature regimes are consistent and indeed the 
diagram in figure 1 is commutative for the considered choice of operators. 

We proceed now similarly and compute the two-point correlation functions in- 
volving various other operators. In what follows we will be less explicit as for 
(J(r)J(O)) in the derivation of the finite temperature and mass correlation function 
and just quote the final results. Thereafter, we carry out the various limits by using 
the formulae quoted in the appendix. 

Recalling that the only non-zero form factors of the trace of the energy-momentum 
tensor 6 are 



) = -2-nim 2 sinh(# - 0)/2, 
we obtain, according to ( |2.13|) , for its mutual correlation function 



(G(r)9(0)) 



T.m 



2m 4 



k£(m, r, T) 2 — Ko(m, r, Tf 



(3.11) 



(3.12) 
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We can verify the commutativity of the diagram in figure 1 similarly as for the 
current- current correlator 



lim 

m— >0 

lim 



L hm<e(r)e(o)) Tim 
Km <e(r)e(o)> Ttm 



lim [2m 4 (Kf(rm) - K$(rm))} = 0, (3.13) 



lim [0] = 



(3.14) 



The conformal limits in (|3.13|) and G3.14Q reflect of course the vanishing of the trace 
of the energy-momentum tensor. Noting that the energy density operator e with 
conformal dimension A e = A 6 = 1/2 is related to the trace as = me, we obtain 
more interesting limits 



lim 

in— >0 

lim 



lim (e(r)e(O)) 



T,m 



lim [2m 2 (K\{r-m) — K^rm))] 



lim (e(r)e(O)) 



.m->0 



T,m 



lim 



2vr 2 T 2 



sin (irrT) 



2A e +2A £ ' 



r 



r, (3.15) 
(3.16) 



which are again consistent. Recalling that for the (++)-component of the energy 
momentum tensor T ++ = T with Af = 2, Af = the only non-vanishing form 
factors are 



? T\ii, 



ni/2m 2 exp(# + 0) sinh(fl - 0)/2. 



(3.17) 



The mutual correlation function results to 



(f(r)T(O)) 



m 



K^(m, r, T)K+(m, r, T) — K 2 (m, r, T) z 



(3.18) 



Once again the commutativity of the diagram in figure 1 is confirmed 



lim lim (T(r)T(0))„ 



lim 



lim<f(r)f(0)) Tm 
We also compute 



lim 



lim 



— ( < K 1 (rm)K 3 (rm) — fT 2 (rm)) 



tt 4 T 4 



2 sin (nrT) 



(3.20) 



f(r)e(0)> T 



m 

"' " ~2~ 



K 1 (m, r, T) 2 — K^m, r, T)i^ 2 h (m, r, T) 



together with the expected limiting behaviour 



(3.21) 



lim 

m— >0 

lim 



hm<f(r)e(0)> 



.T-t-O 



lim 

m^O 



m 



(fT 2 (rm) — i^"o( rm )-^2(^^T.)) 



lim <r(r)9(0)> 

m— >0 



lim fOl = . 



0,(3.22) 
(3.23) 
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Once again, the conformal limits in ( |3.22| ) and ( |3.23| ) reflect the vanishing of the 
trace of the energy-momentum tensor. Replacing — > me will only change in (|3.22|) 



m 



and the remaining limits are the expected ones. 



However, there are also operators for which the prescription does nor work so 
smoothly. As an example, we now want to compute (T(r)/x(0)) Tm , with /i be- 
ing the disorder field of the complex free Fermion theory. For this purpose, it is 
necessary first to recall the expressions of the form factors related to this field, 
which were computed in [p3|, and shown to be different from the ones associated to 



the counterpart of this field in the Ising model. Similarly as for the latter model, 
it was found that only the form factors involving an even number of particles are 
non-zero. Since for T the only non-vanishing form factors are the two-particle ones, 
the only non-vanishing form factors of the field [i which will be relevant for these 
computations are 



-9, -6) = i/2 (/i) T=0 exp[(0 - 0)/2] cosh" 



With these data we compute, again according to ( |2.13| ) 
<f(r)/i(0)) r , 



l/T 

T 3 (M)r=0 f dte -2tm 



16 



2m / ml . 
— $(-e",_ tT)<L> 



0)/2 . (3.24) 



-e-r,-,tT 



-<1> j - e --,-,tTj +3$ f -e--,-,tTj $ I - e ~-,-,tT ) / 



1 



. (3.25) 



We employed Lerch's transcendental function <f>(x,s,ct) (see appendix). In com- 
parison with our previous computations we have one integration remaining. This 
results from the fact that unlike before we have now a term (cosh# + cosh^) in the 
denominator, which we eliminate by a differentiation with respect to r and a subse- 
quent integration. Alternatively, we obtain the same result by a direct computation 
using a variable substitution similar as in [20| sinh#/2 = rcos0, sinh#/2 
. By employing ( |A.14j ) the limits come out as 



r sin i 



lim 

m— >0 



hm(f(r)M0)) 



,T-»0 



T,m 



lim 

m— >0 



<A*> 



-2rm 



T=0 



16r 2 



<A*> 



T=0 



16r 2 



(3.26) 



However, starting with the limit m — > in ( |3.25|) is problematic, because the 
expression lim^-i $(x, s, a) for Res < 1 is not well defined. We obtain similar 
phenomena for other correlation functions involving different operators, such as 

In principle one could extend this list of correlation functions involving various 
other operators and support more and more our overall conclusion, namely that the 
conjecture of a dressed form factor expansion fl2.13| ) is meaningful for free theories 
even when the underlying statistics is anyonic. Next we want to see whether the 
picture still remains the same for dynamically interacting theories. 
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4 The scaling Yang-Lee Model 

r<y\ 

The scaling Yang-Lee model (or minimal A 2 -afnne Toda field theory), like its 
conformal counterpart with Virasoro central charge c = —22/5, is one of the simplest 
interacting integrable quantum field theories in 1+1 space-time dimensions. It is an 
ideal starting point to test general ideas, since it is comprised of only one massive 
particle which couples to itself. This is reflected by the pole in the physical sheet 
of its scattering matrix 

S YL(9) = ^ e + i ™' 3 (4.1) 
sinn v — i sin tt/S 



which was proposed in [p5 |. Closed formulae for all n-particle form factors for 



various components of the energy- momentum tensor were computed in P6J. For 



our purposes we will just require the ones up to two particles. We recall from [gC 
in a slightly different notation the form factors associated to T ++ = T 

F ° = -7T' F ^ ) = -^i7i e (4 - 2) 

F f (f) 0^ 7rm 2 coshfl 12 -l . , ( 

where the so-called minimal form factor is 



Fmm(d) = exp 

and v is a constant given by 



oo 



dt sinh | sinh | sinh / t(m — 9) 



sinh 2 1 \ 2tt 



(4.4) 



/ f°° dt sinh | sinh | sinh | \ 

u = exp 2 / 2 ^ " =111154... (4.5) 

V Jo t sinh t J 

Our aim is to use these expressions and compute by means of the dressed form 
factor formula (|2.13|) the two-point correlation functions. Unfortunately, to our 



knowledge there exists no computation in the massive and temperature dependent 
situation to compare with. However, in the massless case we have two benchmarks, 
namely 

(T(r)T(0)) M = -HJ. (4.6) 

( f (r)T(P)) M = -t • ^r-r\ ■ (4 - 7) 

,m_u 5 sm [irrl ) 



Here ( |4.6| ) is just the well known two-point function from conformal field theory 
(T(r)T(0)) T _ m _ = c/2r~ 4 , with c = —22/5, and fl4.7|) is this formula mapped to 
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the cylinder according to ( |2.21| ) with = 2. Let us therefore compute the mass- 
less scattering matrices and form factors from (|4.1| ) and ( |4.2|) , ( [4.3|) , respectively. 
According to the prescription ( 2.24Q outlined at the end of section 2, we compute 



Srr(9) = S LL {9) = S YL (0), S rl (9) = S LR (9) = 1 



and 



if 



2,LL 



r 2,LR 



2,RL 



m m 



.20 



it m 



- cosh 9 



12 



2 cosh# 



12 



1/2 



7 12J 



(4.8) 

(4.9) 
(4.10) 

(4.11) 



These expressions also exemplify our remark at the end of section 2, namely, that 
one can not take the scattering matrices ( [O ) and compute the form factors there- 
after from first principles. In that case we would obtain an L and LL contribution, 
which are evidently vanishing when we carry out the method directly on the formu- 
lae ( |4.2|) and ( |4.3| ) for the massive regime. The reason is simply that the "massless 
prescription" is only a way to carry out the limit starting with the massive expres- 
sions, but not a first principle concept. Nonetheless, viewing it in this sense it works 
extremely well. 

Let us commence with the zero temperature case. The one-particle contribution 
can be computed analytically 



(f(r)f(0)> 



(i) 

T=m=0 



d9 
2^ 



71 



V3 1 



2.2020498 



2v 2 r 4 



(4.12) 

We observe, that the expected value is already almost saturated. To compute the 
two particle contribution is a fairly simple numerical exercise. We evaluate 



(f(r)f(O)) 



(2) 

T=m=0 



d9 x d9 2 
2(2tt) 2 



F2,Rr( 



h,9 2 



-rhr(e e i+e e 2) 



(4.13) 



and present our results in table 1 



r 


0.0001 


0.0003 


0.0005 


0.001 


0.01 


0.1 


r*(T(r)T(0)>^ m=n 10 3 


0.5631 


1.9435 


2.0457 


2.0487 


2.0487 


2.0487 



Table 1: Two-particle contribution to (T(r)T(0)) T= _ Q . 

As mentioned in ||26|| , when summing up according to (|2.13|) with /^(^,T) = 1, 



this contribution enters with a positive sign in comparison with the one-particle 
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contribution due to the non-unitarity of the model. As expected, similar to the 
(T = 0, m ^ 0)-case, carried out in we observe an extremely fast convergence 
of the series towards the expected value (|4.6|) . The numbers in table 1 confirm the 
general observation, which was also made in , that for extremely small values of 
r the higher particle contributions become more important. Considering a regime 
for r > 0.001, it will not be necessary to include also the three-particle contribution 
in order to reach our main conclusion. Nonetheless, in principle this could be done 
easily with some Monte Carlo integration, just at the cost of longer computing time, 



since m 



26| all n-particle form factors were already provided. 



Notice further that the parameter rh plays no role anymore. As observed before, 



in the analytical computation ( 4.12|) it cancels explicitly. This phenomenon is less 
apparent in the two-particle formula ( 1.13 ), but we convinced ourselves that rh may 
be re-scaled without altering the outcome of the numerical computation. Thus the 
expressions are mass independent as they should be. In what follows we will use 
this fact to simplify notations and scale rh to one. 

Let us now embark upon the non-zero temperature case. By means of the 
conjecture (|2.13 ), we have to compute for the one-particle contribution 



;f(r)f(0)> 



(i) 

T,m=0 



de 

+f-{0,T) 



2 



-rTe(6,T) 



e rTe(8,T) ^6 



(4.14) 



According to the conjectures in |7], [Kj or if we extend the proposal in || to the 
two-point function, we can choose for the s(9,T) in ( |4.14j ) and the corresponding 
dressing functions (|2.15 ) either £tba(#, T), determined by the massless version of the 
TBA-equation ( pTTBj ) or Ef Tee (9,T) = e 6 /T, respectively. Solving first numerically 
the massless TBA-equation by means of a standard iteration procedure, we can 
compare the two dressing functions /±(£tba(#, T)) and f±(sf ree (9, T)). Our results 
are depicted in figure 2. 
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Figure 2: Two possible dressing functions /±(£tba(^> T)) and f±(E{ ree (9,T)). 

We plotted the corresponding functions for the left and right movers in order to 
exhibit the symmetry of the TBA solutions. The solutions for the left movers are 
not important for what follows. We observe that f_(e^ B j^(6,T)) acquires plateaux 
(0,T)) = 2/(1 + ^) which are characteristic for all minimal affine Toda 
field theories. Furthermore, we see the well-known fact that in the large rapidity 
regime the TBA solutions merge with the free one. Nonetheless, the two functions 
f± BA (9,T) and f± cc (8,T) differ quite substantially within a large rapidity regime. 
However, in the computation in which they are actually needed, namely in (|4.14 ) 



this regime is negligible. This is essentially due to the factor e 46 \ This means the 
issue of controversy on the difference between the two functions f± BA (9,T) and 
f± ee (8,T) is irrelevant if we would extend it to the context of the two-point func- 
tions. Let us therefore take f± ee (8, T), for which we can compute ( f4.14| ) analytically 



T(r)T(0)>« m=0 = -^T 4 [$(-1, 4, rT) + $(-1, 4, 1 - rT)] . (4.15) 

Here s, a) is again Lerch's transcendental function, which was already encoun- 
tered in section 3 (see also appendix). Extrapolating the behaviour from the cases 
(m^0,T = 0) and (T = 0,m = 0), we expect that the main contribution to the 
sum comes from the one-particle form factors. It is therefore instructive to com- 
pare the ratio of the expected value ( [Op and (T(r)T(O)) * 0m _ Q . We depict this 
comparison in figure 3. 
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Figure 3: Exact correlation function versus dressed one particle form factor contribution 
G(rT):=(f (r)f (0)>^ Ojm=o / (f(r)f (0)>« 0m=0 . 

We observe a deviation of up to 300%, which when recalling the excellent agree- 
ment at this level of the cases (m 7^ 0,T = 0) and (T = 0, m = 0), sheds a 
rather pessimistic light on the working of the conjecture ( |2.13| ) in the interact- 
ing case. In fact, we only observe a reasonable match when < rT < 0.01 or 
0.99 < rT < 1.0, but this is unfortunately just a regime in which we can ap- 
proximate in ( }4.7| ) sin (nrT) ~ (7rrT) 4 such that the temperature effect becomes 
irrelevant. For the sceptical reader we include in figure 3 also some points obtained 



numerically be using f± BA (9,T) instead of f± ee (9,T) in ( |4.14| ). The two different 



cases may hardly be distinguished. 

Let us see whether the next order contributions can improve the situation. For 
this we have to compute 



;f(r)f(0)> 



(2) 

T,m=0 



d9 x d9 2 



2(2tt) 2 
+/_(0 1) T)./L(0 2 ,T) 

+2/ + (0 ls T)/_(0 2 ,T) 



U{9 U T)U{9 2 ,T) 



FT 



7 Irr{9i, ®2 



-rT(e(6uT)+e(0 2 ,T)) 



3 rT(e(0i,T)+ e (0 2 ,T)) 



2,RR 



(01,02 



Z7T 



-rT{e{e 1 ,T)-e{9 2 ,T)) 



.(4.16) 



We have already all the ingredients to compute this apart from the particle-hole 
form factor -F^^(0i, 9 2 + Z7r). We compute 



2,RR\ 



h,e 2 



7T 



2coshg 12 + 1 \ tanh 2 (g 12 /2) c g i+g2 



2z/ 8 \2cosh# 



12 



(4.17) 



7 12 
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Assembling all we find similar values in the entire range of < rT < 1 

(f (r)T (0)>^ 2) m=0 / <f (r)f (0)> w ~ 0.001 ± 4 x 10^ 4 , (4.18) 

independently of the different choices for the dressing functions. Thus, assuming 
that the convergence of the series in ( |2.13|) does not change radically when the 
temperature is switched on, the higher order n-particle contributions will not rescue 
the proposal for this case. 



5 Conclusions 

We provided more evidence which supports the proposal of LeClair, Lesage, Sachdev 
and Saleur || to use dressed form factors for the computation of two-point corre- 
lation functions. The method seems to work well for the free Fermion case and in 
addition for theories with anyonic statistics. Concerning dynamically interacting 
theories we reach a similar conclusion as drawn in [BJ on the base of an example 
involving a chemical potential: Namely that it fails to work. As a simple counter 
example we have studied the scaling Yang-Lee model. This conclusion is reached in- 
dependently of the choices for the dressing functions /(£tba($> T)) or f(ef vee (9, T)). 

Despite this slightly pessimistic result concerning the proposal in its present 
form, it was shown in that it also works for the computation of correlation 
functions in the presence of boundaries, albeit only for the Ising model. Based on 
this result one may conjecture that it can also be successfully applied to defect 



systems p7 |. In fact in this context the only interesting integrable theories are 
those for which the approach seems to work for the bulk theories, namely free 
theories, possibly with anyonic statistical interaction. It was shown recently p8 |, 



that these theories are the only bulk theories which allow a simultaneous occurrence 
of reflection and transmission. 

At present the challenge remains to show whether the form factor approach to 
compute correlation functions can be extended successfully to the finite temperature 
regime in complete generality, namely also for dynamically interacting theories. In 
order to complete this task it would also be interesting to compare with existing 



alternative approaches, e.g., [29 
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Appendix 

In this appendix we assemble some properties of various functions which are impor- 
tant for our computations. Some of them are standard whereas others are specific 
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to the present context. One of the most ubiquitous functions in this context are 
the modified Bessel functions, whose integral representations are given by 



KJz) 



7T 



dt exp(— zcosht) coshat for |argz| < — . 



We recall the well-known limiting behaviour 
lim K a {x) ~ 2 a - 1 T(a)x- a Re a > 0, 

x— >0 

It is convenient to introduce the function 

oo 

^(m,r,T)^(-ir 



lim K (x) ~ — In x . 

x^Q 



n=0 



/nm 

K a [ — + rm) ±K 



\ T 



[n + l)m 



T 



— rm 



(A.1) 



(A.2) 



(A.3) 



which will appear as a building block in the computation of many finite tempera- 
ture correlation functions. Since we intend to investigate the commutativity of the 
diagram in figure 1, various limits of this function will be required frequently. For 
Re a > we compute with ([A.2 ) the massless limit 



Via) x 



lim^(m,r,T)~-gJ>l 

n=0 



(nm \~ a f(n + l)m 

1- rm. + - - — 



+ rm ± — rm 

V T J \ T 



(A.4) 

The sum can be evaluated explicitly. We just report on the cases which are impor- 
tant for our analysis 

7rT 1 

lim K+(m,r,T) — , (A.5) 

m-,o m sm(7rri ) 

T -v / ^ ( \ 2 2 cot (ivrT) 

lim K-(m,r,T) ~ — . 2 / ' , A.6 
m-»-o \m J sin (7rri J 

/vrT\ 3 6 + 2cos(27rrT) . A „, 

lim4(m,r,T) ~ — . 3 / — ; , (A.7) 



m^o \m J sm s (wrT) 

lim m 4 K^(m,r,T) ~ 0, (A.8) 

lim m A K£(m,r,T)K£(m,r,T) ~ 0. (A.9) 

The zero temperature limit is more easily computed, just by noting that in the sum 
of (|A.3|) only the n = term survives 

lim k±(m,r,T) ~ (rm) . (A. 10) 



19 



A further function which frequently occurs is Lerch's transcendental function, whose 
sum representation is 

°^ x n 

$0, s,a)=y 7 r- for Id < 1, a 4. Z«. (A.ll) 

n=0 V ' 

In many cases, however, we require precisely the value x — > — 1 in ( A.ll ). The 



convergence problem can be circumvented by exploiting the fact that in the limit 
x —>■ — 1 we can express <&(x, s, a) in terms of Riemann zeta functions 



C(s,a) = y \ for Res>l, (A.12) 

^— » In -4- n/ I s 



n In + a, 

n=0 v ' 

instead 

lmi i s, a) = i [C(s, a/2) - ((s, (1 + a)/2)] . (A.13) 
This leaves the problem lim^_i s, a) for Res < 1. Also the limit 

lim $(-e- 1/Q , s, a) ~ a" s , (A.14) 

is needed in section 3. 
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